Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



1 



A new part-per-million measurement of the positive muon lifetime and 
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The Fermi Constant, Gf, describes the strength of the weak force and is determined most precisely 
from the mean life of the positive muon, r^. Advances in theory have reduced the theoretical 
uncertainty on as calculated from T/^ to a few tenths of a part per million (ppm). Until recently, 
the remaining uncertainty on was entirely experimental and dominated by the uncertainty on 
T/j,. We report the MuLan collaboration's recent 1.0 ppm measurement of the positive muon lifetime. 
This measurement is over a factor of 15 more precise than any previous measurement, and is the most 
precise particle lifetime ever measured. The experiment used a time- structured low-energy muon 
beam and an array of plastic scintillators read-out by waveform digitizers and a fast data acquisition 
system to record over 2 x 10^^ muon decays. Two different in-vacuum muon-stopping targets were 
used in separate data-taking periods. The results from these two data-taking periods are in excellent 
agreement. The combined results give (MuLan) = 2196980.3(2.2) ps. This measurement of the 
muon lifetime gives the most precise value for the Fermi Constant: (MuLan) = 1.1663788(7) x 
10~^ GeV~^ (0.6 ppm). The lifetime is also used to extract the /j.~p singlet capture rate, which 
determines the proton's weak induced pseudoscalar coupling gp. 

I. INTRODUCTION 

The Standard Model of particle physics has given an excellent description of particle interactions since it was 
first formulated in the early 1970's. The predictive power of the Standard Model depends on its well-measured 
input parameters. These parameters include the masses of the leptons and quarks and the mixing angles of the 
unitary CKM quark mixing matrix [9]. The Standard Model also describes three of the fundamental forces of 
nature: the electromagnetic, weak, and strong forces. 

The weak force requires three parameters to characterize its interactions. The most precisely measured 
parameters are the Fermi Constant Gf, the mass of the Z-boson M^, and the fine structure constant a. The 
Fermi Constant describes the strength of the weak force, and is related to the electroweak gauge coupling g by 

where Ar includes all higher-order weak interaction loops. 

The Fermi Constant is determined most precisely from muon decay. Muon decay is most likely to proceed 
via the reaction /i+ ^'^^e^^- In Fermi theory, muon decay is described by a four-fermion contact interaction. 
The relation between the muon lifetime and the Fermi Constant is 

1 Glml 

where includes phase space, QED, hadronic, and radiative corrections. In 1999, the 2-loop QED corrections 
were calculated, which reduced the theoretical uncertainty on the Fermi Constant as calculated from to less 
than 0.3 ppm; previously it was the dominant uncertainty [1]. This improvement in theory motivated new 
measurements of the muon lifetime: MuLan [2j and FAST 



II. EXPERIMENT OVERVIEW 



The Muon Lifetime Analysis (MuLan) experiment used a time-structured surface- muon beam and a symmetric 
detector to observe over 2 x 10^^ muon decays. The experiment took place in the ttES beamline at the Paul 
Scherrer Institute (PSI) in Villigen, Switzerland. PSI has the most intense continuous proton beam in the world. 



* This work was performed while the author was at the University of lUinois at Urbana- Champaign, USA 
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The > 2 mA, 590 MeV proton beam is directed through a 4-6 cm carbon target, producing pions. The positive 
pions that decay at rest near the surface of the carbon target produce positive muons which are cohected by 
the 7rE3 beamhne and guided toward the MuLan experiment. 

A time structure of 5 /is beam-on and 22 /is beam-off is imposed on the continuous muon beam by a 25 kV 
electrostatic kicker [4]. During the 5 /is beam-on period, muons propagate down the beampipe and stop in a 
thin target at the center of the MuLan detector. During the 22 /is beam-off measurement period, the muons 
are deflected into a beam collimator upstream of the experiment. The beam extinction, defined by the ratio of 
beam-on to beam-off muon rate at the detector, is typically around 1000. 

The muons stop in a target at the center of the MuLan detector. The negative helicity of the muons is 
preserved during beam transport, and muons stop in the target with ~ 100% polarization. A slow precession of 
the ensemble muon spin from the Earth's magnetic field would introduce an oscillation in the observed muon 
decays which could pull the lifetime fit. The muon stopping targets were chosen to dephase or precess the muon 
ensemble to remove or obviate the muon spin precession. During the 2006 run, the target was composed of 
a ferromagnetic alloy of ~ 60% Iron, ^ 30% Chromium, and ~ 10% Cobalt, called Arnokrome^^ III (AK-3) 
[5[. The 0.4-T internal magnetic field precessed the muon spins with period 18 ns. This precession dephased 
the muon ensemble during the 5 /is beam-on period. During the 2007 run, a crystalline quartz target with 
130-Gauss externally-applied magnetic field was used. In quartz, the majority of the muons form a muonium 
bound state consisting of a muon and an electron. The magnetic moment of the paramagnetic bound state is 
^ 200 X higher than that of the muon alone. Muons in the paramagnetic state precess with period 2.6 ns, and 
muons in the diamagnetic state precess with period 550 ns. The 550 ns muon precession period is observed. 

The MuLan detector symmetrically surrounds the muon stopping target and records the outgoing positrons 
from muon decay (Fig. [1]). The detector is composed of 170 triangular scintillator tile pairs. The coincident 
signal between inner and outer tiles of a pair defines a through-going particle. The scintillation light from 
each tile is collected on one edge of each tile and guided toward a 29-mm photomultiplier tube (PMT). The 
analog signal from each PMT is read-out by a 450 MHz waveform-digitizer (WFD). When the signal goes over 
threshold (a "hit"), 24 waveform samples are recorded along with the channel and the time. Every hit is stored 
to tape, and ^ 60 terabytes are required to store the information from 10^^ muon decays. 



FIG. 1: The MuLan detector. Muons are transported in vaccuum to a stopping target in the center of the 
detector. The stopping target was composed of AK-3 in 2006 and crystalline quartz in 2007. An external 
magnetic field was applied to the quartz target using a Halbach array of permanent magnets. Both targets can 
open to allow the beam to pass to a wire chamber downstream of the experiment (EMC). 

In the analysis, coincidences are identified as hits in the inner and outer tiles of a pair within a given time 
interval, usually 11 ns. These coincidences are histogrammed versus time since the beginning of the measurement 
period (Fig. [2]). This lifetime histogram is then fit with function 



where Nq is a normalization constant, is the muon lifetime, 5 is a flat background, S{t) is a small oscillation 
from the readout electronics (see section UlI Ap . and A is the scale of the electronics oscillation. During the 




f{t) = (1 + AS{t))Noe 




(3) 
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analysis the lifetime was represented by defined as a ppm- difference from a reference value: 

= TO (1 + i?/10^) . (4) 




To prevent experimenter bias during the analysis, the WFD digitization frequency was blinded and the muon 
lifetime was known only in "clock ticks" (ct). After the analysis of the systematic uncertainties was complete, 
the clock frequency was unblinded to reveal the measured muon lifetime. 



III. SYSTEMATIC UNCERTAINTIES 



During the 22- jas measurement period (about 10 muon lifetimes), the hit rate in the detector changes by 3 
orders of magnitude. Any rate- or time-dependent change in detector performance during the measurement 
period could distort the fitted muon lifetime. A list of systematic and statistical uncertainties is given in table 
[D Two of these uncertainties, gain stability and pileup, will be discussed in detail. 



TABLE I: Systematic and statistical uncertainties in ppm units. The errors in different rows of the table are 
not correlated to each other. Where only one error appears in a given row, the effect is 100% correlated 



between the two run periods. 


Effect uncertainty in ppm 


R06 ROT 


Kicker stability 


0.20 0.07 


Spin precession / relaxation 


0.10 0.20 


Pileup 


0.20 


Gain stability 


0.25 


Upstream muon stops 


0.10 


Timing stability 


0.12 


Clock calibration 


0.03 


Total systematic 


0.42 0.42 


Statistical uncertainty 


1.14 1.68 



A. Gain Stability 



A change in signal height, or gain, will change the number of hits over threshold. If the gain changes in time- 
dependent way it will introduce a perturbation on the lifetime histogram which may pull the fitted lifetime. 
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For example, Fig. [3] shows a small oscillation in the pulse height most-probable- value early in the measurement 
period, which originates in the electronics. Multiple analysis techniques were employed to characterize the effect 
of this oscillation on the fitted muon lifetime. When the oscillation is taken into account, either by including it 
in the fit function or correcting the lifetime histogram, the fitted muon lifetime changes by 0.50 ± 0.25 ppm. 
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FIG. 3: Most probable value for the pulse height (MPV) vs. time in measurement period. The amplitude 
spectrum of the background, determined in the time range 9500-9600 ct (20.9-21.1 /is), was subtracted. An 
electronics oscillation is visible early in the measurement period. Inset: amplitude minus threshold, summed 

over all detector PMTs. At threshold, the fractional change in counts for a 1 ADC change in threshold is 
^ 3.0 X 10~^. For comparison, at threshold+40 ADC, the fractional change in counts is ~ 114.4 x 10~^. The 
effect on the lifetime histogram is found by scaling the oscillation in MPV by the fractional number of counts 
at the amplitude threshold. Taking the oscillation into account shifts the fitted muon lifetime by 

0.50 ±0.25 ppm. 



B. Pileup 

Although the waveform digitizers (WFDs) provide complete waveforms for each hit, the analysis software 
cannot reliably distinguish two hits separated in time by less than 5 WFD samples (^11 ns). When two hits 
occur within this so-called deadtime (DT), the second hit is lost. This effect is called pileup. 

Pileup is statistically corrected using the data. When a hit is observed at time ti in beam filH, a search 
is triggered in the next beam filH + 1 from time ti to ti + DT. If a hit occurs in the same channel in this 
search window, it is stored in a pileup-correction histogram. The pileup-correction histogram is added to the 
main lifetime histogram prior to fitting. This "shadow-window" method statistically corrects for the pileup 
by counting certain hits twice to compensate for the hits which are lost. Several higher-order pileup effects 
were also considered, including triple-hits, accidental coincidences, and "softening" of the coincidence window 
boundaries due to jitter in the inner-outer coincidence times (Fig. |4]). 

The full pileup correction method was tested in Monte-Carlo with 10^^ simulated muon decays (Fig. [5a|). The 
correction shows excellent consistency with the simulated lifetime over a range of artificially applied deadtimes 
(ADTs). When this technique is applied to the data (Fig. [5b]), a residual slope exists of order 0.008 ppm/ns 
of applied deadtime, corresponding to a ~ 0.1% pileup undercorrection. This undercorrection is attributed to 
fluctuations in the PSI proton beam, resulting in small variations in muon rate from one beam fill to the next. 
The extrapolation to ADT=0 gives the true muon lifetime. 
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FIG. 4: Pileup correction histograms. Even negligibly small O(10~^^) effects were considered. Note that 

lifetime fits are approximate. 




FIG. 5: a) Monte-Carlo pileup reconstruction for 10^^ simulated muon decays. The simulated lifetime R value 
is shown at the left for ADT=0 with statistical error bars 1.19 ppm, larger than the range of the plot. The 
other points at selected ADTs are shown with error bars allowed by the deadtime correction. The difference 
between the simulated lifetime and the average reconstructed value is 0.16 ppm. b) R vs. ADT for the 2006 

dataset, after all deadtime corrections have been applied. The error bars are the allowed deviations in R from 
the statistical deadtime reconstruction, and are smaller than the overall statistical uncertainty of 1.14 ppm. 



IV. RESULTS 



In addition to the lifetime fit described by equation [3] and shown in Fig. [21 an exhaustive set of cross-checks 
was performed. One of the most powerful cross-checks is a plot of the fitted muon lifetime vs. the fit start time. 
As the start time of the fit is increased, any potential early-time effects are excluded from the fit. A trend in 
the plot of lifetime vs. fit start time is an indication of an undiscovered systematic effect. This plot is shown 
in Fig. [6] for the 2006 and 2007 run periods. The error bars are statistical, and the bands show the allowed 
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1-sigma statistical variation in the fit. No systematic trends are seen. 
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FIG. 6: Consistency of lifetime fit vs. fit start time. Statistical errors are shown for each point, and the bands 
show the 1-sigma allowed variation. No systematic trends are seen in either the 2006 or 2007 run periods. 

During data analysis, the muon lifetime was known only in clock-tick units, and the blinded clock frequency 
was different for the two running periods. After all systematic cross-checks were complete, the electronics clock 
frequency was unblinded to reveal the measured muon lifetimes: 

r^(R06) = 2196979.9 ± 2.5 ± 0.9 ps, 

r^(R07) = 2196981.2 ± 3.7 ± 0.9 ps. (5) 

where the first uncertainties are statistical and the second are systematic. The results from the two running 
periods are in excellent agreement. Combining these results yields 

r^(MuLan) = 2196980.3 ± 2.2 ps (1.0 ppm). (6) 

These results are shown in the context of previous precision muon lifetime experiments in Fig. [71 
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FIG. 7: Precision muon lifetime measurements. Results from the 2006 and 2007 MuLan run periods are 
plotted separately to show their consistency. The gray band shows the combined MuLan result [7,] . 

From this value for , we derive a new value for G p following equation O 

GF(MuLan) = 1.1663788(7) x 10"^ GeV"^ (0.6 ppm), (7) 

where the dominant contribution to the Gp uncertainty is 0.5 ppm from the muon lifetime. 

Aside from a new determination of Gf-, is also used in muon capture experiments. The MuCap experiment 
measures the lifetime of negative muons in protium gas [6]. The difference between the positive and negative 
muon lifetimes is used to extract the singlet capture rate for the process jii~p u^n. The capture rate is then 
used to derive the pseudoscalar form factor of the proton, gp. Fig. [8] shows that using the MuLan result for 
shifts the value for gp into better agreement with theory. 

V. SUMMARY 

The Muon Lifetime Analysis (MuLan) experiment has determined a new precise value for the positive muon 
lifetime, r+. The experiment measured over 10^^ muon decays using a time-structured surface- muon beam and 
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FIG. 8: Results from the MuCap experiment, which measures the negative muon capture rate on the proton 
and derives the pseudoscalar form factor of the proton, gp. A precise value of the positive muon lifetime is 
necessary to derive the capture rate from the measured negative muon lifetime in hydrogen. Using the MuLan 
result shifts the derived value of gp into better agreement with theory. The parameter Aop is the transition 
rate of muonic molecular hydrogen (p/ip) from the ortho to the para state. These ortho and para states have 
different muon capture rates. MuCap is designed to be relatively insensitive to the transition rate, which 
allows a clear extraction of gp from the measured muon capture rate. 



a symmetric detector. The result, r^(MuLan) = 2196980.3 ± 2.2 ps (1.0 ppm), is more than 15 times more 
precise than any previous measurement and is used to provide the most precise value for the Fermi Constant, 
(MuLan) = 1.1663788(7) x 10~^ GeV~^ (0.6 ppm) [8]. The muon lifetime is also used to extract the muon 
capture rate v^n and the pseudoscalar form factor of the proton, gp. 
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